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ABSTRACT 


The major driver for broadband wireless communications has been reliable, high-data 
rate services. In wireless communication, the multipath fading constitutes a bottleneck 
for increasing data rates and causes performance degradation. To combat fading, we can 
use diversity. Wireless systems with multiple antennas at the transmitter and receiver 


have much larger capacity in fading channels than standard wireless systems. 


The objective of this thesis is to investigate the transmission scheme provided by 


matrix A and B in the 802.16 standard and show how it can be implemented. 


The research focuses on using maximal-ratio combining (MRC) to demodulate 
the transmitted symbols. Modifications to the existing matrix by using more frequency 
bands were introduced; this reduces the number of transmitting antennas and uses fewer 
time slots to transmit the same number of symbols. The modulator and demodulator 
design is also discussed. The performance of orthogonal and non-orthogonal space time 


codes (STC) are evaluated. 
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EXECUTIVE SUMMARY 


Wireless systems with multiple antennas at the transmitter and receiver have a much 


larger capacity to combat fading channels than standard wireless systems. 


These systems use diversity techniques to improve the reliability of transmission 
and reduce the effect of multipath fading without increasing the transmitted power or 
increasing the bandwidth. Unlike time and frequency diversity, space diversity does not 
require more bandwidth. This property makes space diversity very popular for high data 
rate wireless communications. This research examines the maximal-ratio combining 
(MRC) method. This method maximizes the signal-to-noise ratio (SNR) since the 
maximum output SNR is equal to the sum of the instantaneous SNRs of the individual 


signals. 


Two matrices, matrix A and B in the 802.16 standard, are evaluated in this 
research. Matrix A refers to the Space time block coding (STBC) technique and multiple 
input multiple-output (MIMO), whereas matrix B refers to the spatial multiplexing 
(SMX) technique. In matrix A, four symbols are transmitted using four antennas in four 
time slots. It has a code rate of one. On the other hand, matrix B transmits eight symbols 
using four antennas in four time slots. Matrix B has a code rate of two. The designs of the 
modulator and demodulator for both matrices are discussed, and the performances are 


evaluated. 


Matrix A is quasi-orthogonal, and matrix B is non-orthogonal. The closed form 
solution to evaluate the performance of matrix A was used in this study, and MATLAB 
was used to evaluate the performance of matrix B. Maximum likelihood (ML) decoding 
was employed to perform the simulations. The simulation results show that there is a 4 
dB drop in performance for quadrature phase-shift keying (QPSK) and 16-quadrature 
amplitude modulation (16QAM) when they are compared to the closed form solution. In 
addition, the comparison of using a single receiving antenna and two receiving antennas 
was performed, and the results show that there is a 11 dB drop in performance when one 


receiving antenna is used. 


xiii 


This research also examined the twelve variations of matrix B and evaluated their 
individual performances. Since matrix B sends out four symbols for the first two time 
slots and the next four symbols in the next two time slots, we can separate the matrix into 
two matrices. Thus, the simulation in MATLAB is a program based on the first two time 
slots and this saves time in the simulation process. The performances for four out of the 
twelve matrices were found to be worse than the other eight matrices. Therefore, we must 


avoid using these four matrices in the transmission process. 


Graph of Bit Error Rate vs E,/N, for QPSK 


2] —:—-- QPSK Closed Form 
| +, 


a Hy 2 
—s—H,., 
—+—H,, 


Bit Error Rate 





E, iN, (dB) 


Figure 1. Simulation results for H)-1, 41-2, H\-3, H\.4 and QPSK closed form solution. 


From Figure 1, the results for H,.4, which is the fourth variation of matrix B, has a 


poorer performance than the rest of the three matrices by 5 dB. 


Finally, we can increase the code rate and have better performance in terms of bit 
error rate (BER) by introducing more frequency bands to the STBC, matrix B. This can 
be proven mathematically using MRC. For example, by introducing four frequency bands 
to matrix B, eight symbols are sent using two time slots, and the code rate is increased 


from two to four. 
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I. INTRODUCTION 


A. OVERVIEW 


The major driver for broadband wireless communications has been reliable, high- 
data-rate services. The scarcity of bandwidth resources motivates researchers to develop 
highly efficient modulation schemes and more efficient coding that utilize less bandwidth 
and, at the same, time provide better quality. 

In wireless communication, multipath fading constitutes a bottleneck for 
increasing data rates and causes performance degradation. Combating this fading requires 
the use of diversity. 

Spatial multiplexing (SMX) using space-time coding (STC) is effective for 
combating fading and enhancing data rates. It relies on simultaneous coding across space 
and time to gain diversity without using too much bandwidth. Transmit diversity has been 
studied extensively as a method to combat impairments in wireless fading channels [1]. 
The main advantage of using transmit diversity is that it is relative simple to implement 
and is feasible to install multiple antennas at the base station (BS). 

One attractive approach of space-time code design is to construct space-time 
block codes (STBCs) from orthogonal designs as proposed by Alamouti [2]. These codes 
achieve full diversity and have fast maximum-likelihood (ML) decoding at the receiver 

The use of maximal-ratio combining (MRC) as a combining method is popular as 
it maximizes the signal-to-noise ratio (SNR) since the maximum output SNR is equal to 


the sum of the instantaneous SNRs of the individual signals [3]. 


B. THESIS OBJECTIVE 


The target of this research is to implement and evaluate the performance of 
WiMAX (Worldwide Interoperability for Microwave Access) orthogonal space time code 
(STC) for orthogonal frequency-division multiple access (OFDMA). The key objective is 
to study the transmission scheme, matrix A and B in the 802.16 standard, and try to 


improve the performances by modifying these matrices. The performances of the STCs 


were evaluated, and their modulator and demodulator designs were discussed. 
Simulations were developed in MATLAB and used for comparison for the respective 


modification of the matrices. 


C; THESIS SYNOPSIS 


The structure of this thesis comprises an Introduction (Chapter I), Background 
(Chapter II), Implementation and Performance Evaluation of WiMAX orthogonal STC 
for OFDMA, an Implementation and Performance Evaluation of WiMAX non-orthogonal 
STC for OFDMA and a Conclusion. The evaluation of matrix A from the 802.16 standard 
and its variations are provided in Chapter II. In Chapter IV, matrix B was evaluated and 
the performance of its variations are presented. Simulations were developed in 
MATLAB. The performance of using one and two receive antennas for decoding was 
compared in this chapter as well. Finally, the conclusion of this thesis report can be found 


in Chapter V. 


Il. BACKGROUND 


The basic principles and background knowledge of diversity [3], the diversity 
combining method [3], SMX and space-time block coding (STBC) are presented in this 
chapter. 


A. DIVERSITY 


Wireless mobile communications use diversity techniques to improve the 
reliability of transmission and reduce the effect of multipath fading without increasing the 
transmitted power or increasing the bandwidth. For example, a time diversity technique 
requires the transmission of same symbols numerous times. Each transmission carries the 


same information but has a small correlation in fading statistics. 


The principle behind diversity is that when two or more independent samples of a 
signal are received, the samples may fade in an uncorrelated manner. Some samples are 
more attenuated, while others are less faded. Therefore, the probability of all samples 
falling below a given level is much lower than the probability of any individual sample 
falling below that level. As such, with proper combinations of the various samples, one 


can greatly reduce the severity of fading and improve the communication link. 


Wireless communication systems use and combine a number of diversity methods 
to achieve the desired performance. Diversity techniques can be classified into time, 


frequency and space diversity. 


1. Time Diversity 


One can employ time diversity by transmitting the same symbol in multiple time 
slots at different times so the receiver will receive uncorrelated fading signals. The 
required time separation is at least the coherence time of the channel, or the reciprocal of 
the fading rate //f7 = c/of-, where fy is the Doppler shift of the received signal, f- is the 


frequency of the carrier, v is the vehicle speed, and c is the speed of light. 


The coherence time is a statistical measure of the period of the time over which 
the channel fading process is correlated. One of the drawbacks of this scheme is that the 


redundancy introduced in the time domain leads to a loss in bandwidth efficiency. 


Zs Frequency Diversity 


Using a number of different frequencies to transmit a message is called frequency 
diversity. For this to work, it is necessary to separate the frequencies enough to ensure 
that each carrier frequency experiences independent fading. To achieve that, the 
frequency separation must be in the order of several times the channel coherence 
bandwidth as this will guarantee that the fading statistics for the various frequencies are 
uncorrelated. The coherence bandwidth is different for different environments. Similar to 
time diversity, frequency diversity will occupy a larger bandwidth as redundancy is 


introduced in the frequency domain. 


3: Space Diversity 


Space diversity, sometimes also known as antenna diversity, uses multiple 
antennas or antenna arrays arranged together in space for transmission and reception. The 
antennas are separated physically by proper distances, so that individual signals are 
uncorrelated. The separation requirements vary with the environment, height of antenna 
and frequency. In space diversity, the repeated copies of the transmitted signals are 
usually given to the receiver in the space domain. Unlike time and frequency diversity, 
space diversity does not require more bandwidth. This property makes space diversity 


very popular for high data rate wireless communications. 


The use of space diversity to implement the STC will be explored in this thesis; 


therefore, a more in-depth discussion of space diversity is necessary. 


One can classify space diversity into two categories: transmit and receive 
diversity. For transmit diversity, multiple antennas are installed at the transmitter station. 
The transmitter processes and combines the messages and delivers them across multiple 
antennas for transmission to one or more receiving antennas. For receive diversity, 


multiple antennas are installed at the receiver station to pick up copies of the transmit 


4 


signals. The replicas of the transmitted signals are properly combined to minimize fading 


due to multipath with the hope to improve the signal-to-noise ratio (SNR). 


Space diversity using different transmit antennas alone cannot provide separate 
copies of the symbol at the receiver since the receive antenna receives the sum of these 
channel-tap weighted copies but does not have the capability to separate them for 
combining. In general, two or more diversity schemes are usually necessary to provide 
multidimensional diversity to meet system performance specifications. A combination of 
two diversity techniques should provide the receiver sufficient statistics to separate 
copies of a transmitted symbol for combining. This is dictated by maximal-ratio 
combining (MRC). As for the receiving end, one can assume that the transmit antennas 
have sufficient inter-element spacing to ensure that all paths from the transmitter to the 


receiver are uncorrelated. 


B. DIVERSITY COMBINING METHODS 


The key feature of all diversity techniques is the low probability of deep fades 
occurring simultaneously in the various diversity subchannels. The performance of the 
communication link depends largely on how the repeated signals are combined, which 
results in an increase to the SNR of the communication system. The type of combining 
methods implemented at the receiver can, therefore, classify diversity schemes. 
According to the implementation complexity and the level of channel state information 
required by the combining method at the receiver, we find that there are four main types 
of combining methods. They are selection combining, switched combining, maximal- 


ratio combining (MRC), and equal gain combining. 


I, Selection Combining 


Consider a receive diversity system with ng receive antennas. The block diagram 
[3] of the selection combining scheme is illustrated in Figure 1. In this system, the signal 
with the largest instantaneous SNR at every symbol interval is selected as the output; the 


output SNR is equal to that of the best incoming signal. In practice, the signal with the 


highest sum of the signal and noise power is usually the one selected, as it is difficult to 
measure the signal and noise power. This method is not recommended when the noise 


power is high. 


Rx | Rx 2 Rx nr 
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Front End Front End 


Logic Selector 


output 


Front End 











Figure 1. | Selection combining method. 


2. Switched Combining 


In switched combining, the receiver scans all the diversity branches and selects 
the particular branch with the SNR above a specified threshold. The receiver then selects 
this signal as the output. When the SNR drops below the threshold, the receiver starts 
scanning again and switches to another branch with the SNR above the threshold. This 


scheme is also called scanning diversity [3] (shown in Figure 2). 


Switched diversity is inferior to selection combining since it does not continually 
pick up the best instantaneous signal. However, it is simpler to design because it does not 


require continuous and simultaneous monitoring of all the diversity branches. 


Both switched and selection diversity schemes obtain their output signals from 
only one of the diversity branches and do not require any prior knowledge of channel 


state information. 
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Figure 2. Switched combining method. 


3. Maximal-Ratio Combining (MRC) 


Maximal-ratio combining is a linear combining method. The received signals are 
individually weighted and added together coherently to obtain a decision statistic. There 
are several ways to choose the weighting factors. The decision statistic is a linear 
combination of a weighted replica of all the received signals. The block diagram [3] for 
maximal-ratio combining diversity is illustrated in Figure 3. In MRC, the weighting 
factor of each receive antenna is chosen to be in proportion to its own SNR. In this 
scheme, each individual signal must be co-phased, weighted with its corresponding 
amplitude, and then summed. This scheme requires knowledge of signal phases and 
channel fading amplitude. This method is also known as optimum combining, since it 
maximizes the SNR. The maximum output SNR is equal to the sum of the instantaneous 


SNRs of the individual signals. 





Figure 3. Maximal-ratio combining method. 


4. Equal Gain Combining 
Equal gain combining is another linear combining method but is suboptimal. It 
does not require the estimation of the fading amplitude for each branch. The receiver 


simply sets the amplitudes of the weighting factors to be unity. 


C. SPACE-TIME BLOCK CODING (STBC) 


Space-time block coding is a technique used in wireless communications to 
transmit multiple copies of a symbol across a number of antennas and combine various 
received versions of the symbol to improve the reliability of the communications link. 


Scattering, reflection, refraction and thermal noise normally affect the transmitted signal. 
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Some of the received copies of the symbol will be ‘better’ than others. This redundancy 
allows for the use of multiple copies to compensate for fading and helps to correctly 
decode the received signal. One can use STC to combine all copies of the same symbol in 


an optimal way to get as much information as possible from each of them. 


iL, Use of Space Time Coding (STC) 


Using an example from [4], we have a space-time block code that allows the 
transmission of m symbol copies in each word of L, symbols; hence, the code rate is m/L,. 
The diversity order of the space-time code is Mz, the number of transmit antennas. To 
separate the symbol copies for combining, it is necessary that the columns of the code 
matrix G be (complex) orthogonal since the transmitter sends all ZL; code words 
concurrently over L, symbol times. One can assume that the channel remains constant for 
m symbol times, and there is perfect channel estimation. Given the transmitted space-time 
block code matrix G, the sufficient statistics for maximal-ratio combining is dictated to 


obtain a matrix H from the vector: 





Gh > Hs 
. 1 
HH =|Al'T, _ 
4T 
where /A =|h, h, .. h,| is the channel tap vector representing L, paths, 
SS [By he ee 1 is the symbol vector representing m transmitted symbols, and H is 


the channel tap matrix. The columns of H must be orthogonal to be able to detect s and to 


decompose s into m separate complex variables. This implies that the columns of G must 


also be orthogonal such that G°G = Is I,,. The matrix H is known by the receiver via 


channel estimation and the known structure of G. 


D. ORTHOGONALITY OF THE CODE MATRIX 


Some STBCs used for WiMAX are orthogonal. This means that the STBC is 
designed such that the vectors representing any pair of columns taken from the coding 


matrix are orthogonal. This enables the design of a simple and optimal decoder at the 


receiver. Most importantly, all but one of the codes that satisfy this criterion must 
sacrifice some proportion of their data rate (the Alamouti code [2] is the only STBC of 
unity rate for complex symbols). 

In addition, there exist quasi-orthogonal STBCs that achieve higher data rates at 
the cost of intersymbol interference (ISI). Thus, their bit error rate (BER) performance is 
lower bounded by the Alamouti unity code that provides ISI free transmissions due to 


orthogonality. 


1. Quasi-Orthogonal STBCs 


The following code is an example of a quasi-orthogonal code: 


Ce. ee, Ce ey 


* * * * 
—C C, =E C. 
2 1 4 ©3 
2 el * * * (2) 
OC, Cy 2 
Cy —C; Cy G 


The code exhibits partial orthogonality, as illustrated in (2). The orthogonality 
criterion only holds for columns (1 and 2), (1 and 3), (2 and 4), and (3 and 4). For 
example, the first two rows of column one and column two will give us a two by two 
matrix. Let us define the two by two matrix as matrix p. Matrix p must have the property 
of (p) p = p(p) =1,, where p’ represents the conjugate and transpose of vector p. The 
same property must apply for row three and row four for them to be considered as quasi- 
orthogonal. All the column pairs mentioned have this property. The code is full-rate and 


requires linear processing at the receiver. 
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Hl. IMPLEMENTATION AND PERFORMANCE EVALUATION 
OF WIMAX ORTHOGONAL STC FOR OFDMA 


First, one should look at matrix A defined in the 802.16 standard [5] to evaluate 
the STC and its variations. One can assume there are multiple antennas on the base 
station (BS) side and one reception antenna on the subscriber station (SS) side. This 
scheme requires multiple-input single-output estimation. Maximal-ratio combining 


(MRC) is used for decoding. 
A. MATRIX A IN 802.16 


1. Matrix A 


Consider the STBC defined by matrix A and given in the 802.16 standard for 
OFDMA (shown in Appendix E): 


Tl T2 T3 T4 
Sy -S, 0 0 ——__ Antennal 
A = Si 5; 0 0 —>___ Antenna2 
7 0 0 S ae Antenna 3 (3) 
3 S4 <= 
* 
0 0 S4 S3 > Antenna 4. 


From the matrix shown in (3), notice that in time slot 1 (T1), symbols s,and s, 
are sent through Antenna | and Antenna 2, respectively. In time slot 2 (T2), symbols —s, 
and s, are sent through Antenna 1 and Antenna 2, respectively. In time slot 3 (T3), 
symbols s,and s, are sent through Antenna 3 and Antenna 4, respectively, and symbols 
—s, and s, are sent through Antenna 3 and Antenna 4 in time slot 4 (T4), respectively. 


This matrix uses four time slots and four antennas to transmit four symbols. Matrix A has 
a diversity of two since it transmits the same symbol twice by using two antennas. For 


example, symbol s, is transmitted by Antenna | at time slot 1 (T1) and Antenna 2 at time 


slot 2 (T2). Matrix A also has a code rate of one since it transmits four symbols in four 
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time slots. The notation s represents the conjugate and transpose of vector s. The 


transmission process at different time slots is illustrated in Figure 4. 


Antl 


oe eee 


fee y 
RF RF 

pe Time Slot 1 Anes 
ae y 
RF RF 





Time Slot 2 








Time Slot 4 


Time Slot 3 Mi 
RF 


Figure 4. Transmission process for matrix A. 


We use MRC to decode matrix A. The first step is to prove if matrix A‘ is 


orthogonal. Matrix A’ is orthogonal only if (A')’4' = 4A’ =J, where A' represents 
the transpose of matrix A, A represents the conjugate of A, and (A')’ represents the 


conjugate transpose of 4" , which is the same as A: 
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(yea SS. 8 : 0 
ss Sy 
S83 
Therefore, 
5,8, + S58, 0 0 0 
ied toes iat 0 5,8, + 558, 0 0 
0 0 S38; + S484 0 
0 0 0 S485 + S48, 
and 
5,8, + S58, 0 0 
Mp age 0 5,8, + S58, 0 
0 0 S353 1 S484 
0 0 0 S48, +5484 
Is," +|s,) 0 0 0 
2 2 
; 0 S| +s 0 0 
AA -A4'A = | H | | ; ; 
0 0 Iss) +|s4| 0 
0 0 0 ssf +|s.f 
Let ||s,|/ =|s,! +|s,/° and ||s, | =|,’ +|s,|° and Q= 4'4" = ATA. 














(4) 


(5) 


(6) 


(7) 


Therefore, 


2: 
























































s, 0 oO O 
0 slr Oo 0 
QO=AA'=A'A = ; 
0 0 |s)? 0 o) 
0 0 0° |s,\P 
a 0 0 0 
Sa 
0 = 0 0 
aT a 
a en (9) 
Is, If 
sal 
ls. 


From Equations (5) and (6), we can prove that matrix A is not orthogonal, but it is 


quasi-orthogonal. 
T 
Next, we can derive the MRC equation to demodulate S = [s, S, 83 s4 
directly from the four antennas. This requires the performance of MRC on matrix A, with 


T 
the receive matrix [Z] =[G][/] + [n], where n= [ 7, coat n,, | is the Gaussian noise 


vector. The A’ matrix is the G matrix for this case. Therefore, 
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Ss; Ss, O 
ge -s, s, 0 
Ol Qe a8 
0 0 -s, 
Sh, + s,h, 
a —sjh, + sh, 
S,h,+ 8,h, 


—S,h, + 3h, 





> > 
3 


nN 

Bl 

h, Nn, 
| n 


(10) 


Next, we make use of the channel estimation to estimate the channel tap vector h. 


Consider an example from [4]. To estimate the channel tap vector h, consider the 


received vector Z=Gh+n where G= A’. Since G’G = Q, the channel estimation can be 


carried out with a known pilot symbol vector s, and G is known to the receiver. The 


sufficient statistics provide the minimum mean-square error (MMSE) of the vector A: 


h=O'Oh+O'G'n 
=h+n. 


(11) 


The jointly complex Gaussian noise vector (Go" )n has the diagonal covariance 


matrix 


E(Q'A" nn’ A'Q") = Q'AE(n’) AQ" 
































=20°0'AA'O" 
=20°Q'QQ" 
=a 0* 
BS: 0 oO 0 
Sq 
0 = 0 0 
5. es 
1 
Is) (12) 
1 
Is) 
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This implies that the components of this noise vector are white and have identical 
variance 2a°Q7'. The components of / can be estimated in the mean square error sense. 
At high SNR, / can serve as the estimate of h. In practice, h may be made available from 


pilot symbols extracted from the pilot tone transmitting simultaneously with the signal. 


Since 
Z, Sh, + sh, n 
Z -sh+sh n 
pee Ne gh + 8, M2 Bie (13) 
Ls S,h,+ s,h, Ny 
Z,| | —syhy + 55h, ny 


the vector Z=Gh+n can be mapped into Y=Hs+N by changing the second element of 
Z=[Z,Z,Z,Z,]' which is Z, = h,s*—hsi+n, into its complex conjugate 
Z; = hjs,—h;s,+n,, and the fourth element, which is from Z, = h,s,—/,s,+n,, into 


its complex conjugate Z; = hs, —h,s,+n, . Thus, 


Zi Sh, + sh, n, 
a ee Z; eh —s,h, +s,h, ; n, 
Zs S,h,+ s,h, n, 
Z; —s,h, +s,h, N, (14) 
h h O O |ls, n 
h, -h 0 0 ; 
=|? 41"? | = CH ILs] + [0] 
0 O h, hy IIs, n, 
0 0 hk, -AIls,| [ni 


From Equation (14), we can generate four separate equations to solve for the four 


unknown symbols, s,,s,,s, and s,. Since the noise 1, m2, 3 and n4 are independent and 
hy, hz, h3 and h4 are known from the channel estimation, we can solve for s,,s,,s,, and 
Si 

The next step is to use the estimated channel taps information [/1, ho, h3, h4] via 
the channel estimation process to estimate /1, h2, h3 and h4, The receiver can then use [H] 


to perform MRC. 
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For simplicity, assume perfect channel estimation where [/, ho, h3, ha] >[h1, ho, 


hz, ha] and apply MRC on Y to obtain the corresponding decision statistics vector LX], 





























































































































where 
2 h 
h, 
X, h; 
XxX h, 
[X]=|  ]= 

Ay é 

x 4 
0 
ae Fy 
h, h, 
Ay ah 
h, h, 
0 0 
0 0 
fy fy 
h,|| \Pa 
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h,|| | 
0 0 
0 0 
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0 |v) 0 0 fs, h, 
= + 
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where ||/,|| and ||/|| are defined as 


[rah = lal + Mal Wal = VMs + [Af (16) 


Next, we would like to determine if the channel is independent. Since 


hh, + hyh, 0 0 0 
was 0 hh, +hyh, 0 0 

0 0 hh, +h,h, 0 
0 0 0 h,h, + h,h, 

la, +\A,/° 0 0 0 (17) 
0 la, +\A,/° 0 0 
0 0 lagf + ry] 0 
0 0 0 lal +|A,) 


letting R = H’H = HH’. Now, 


2 


























R= HH =H'H = ; 18) 


oS 
So 
= 
oS 



































and 
- ; 0 0 0 
1 
0 eer 0 0 
ROS I 
0 1 0 (19) 
al 

0 0 0 L 5 

|| 


From (19), we see that the channel is independent. 
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The block diagrams [5] for the transmitters and receiver are shown in Figure 5 and 6, 


respectively. 
Subchannel IFFT > 
Modulation Input ay IFFT Filter DAC 
Packing nm 
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IFFT Filter DAC 
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Figure 5. | Modulator circuit from the base station for matrix A. 






Diversity Sub- Log- decoder 
Combiner channel likelihood 
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Figure 6. | Demodulator circuit at the subscriber station for matrix A. 


There are four transmitting antennas on the BS side and one receiving antenna on 
the SS side. Two different antennas transmit two different data symbols at the same time. 
The receiver design will have a higher complexity compared to the transmitter. The 
transmitter can be less complex if the matrix can be split into two Alamouti [2] codes. 
This will be discussed in the next design. 

Since both antennas transmit at the same time, they share the same local 
oscillator. Thus, the received signal has exactly the same autocorrelation properties as it 
would for a single antenna. This scheme requires multiple-input single-output (MISO) 


channel estimation. 


20 


From matrix A, we can further break down the matrix into two independent 
; 5, —S, 5, —S, . 
matrices G, = , {and G, = . |that follow the Alamouti STBC. 
Sy S| Ba 83 
From [2], the Alamouti STBC is specified by matrix G; whose codewords 


Ee -s; | and [ 8, 5, ] are the rows of G;. Taking G; as an example, we have 


Antl Ant 2 
i (20) 


T 5S; | S, |—>Symbol Time 1 
G, = 


—s, s, |—>Symbol Time 2 


The columns of G,'are orthogonal and so are the rows; for example, 


G,G,=G,G; = |s| I, where s=[s, s,]'. The Alamouti code allows the transmission of 


two symbols over two antennas in two consecutive symbol times. The code is perfect 


STBC and requires no bandwidth expansion. The sum of the channel-tap weighted code 


eal Sle a 


=| ee Lad | (21) 
hs, —hs, ny 


=Gh+n. 


words plus noise is 


The vector Z=Gh+ncan be mapped into Y= Hs+nby changing the second 
element of Z=[Z, Z,], which is Z,=h,s/—hs,+n,, into its complex 


conjugate Z> = h;s, —h,s, +n,. Thus, we have 


hie Z 7 fas, i N, 
Z, hs, —hs, N, 
= h, h, || 5, = N, 
hy —h, Sy N, 


where N =[N, N,]' =|n, n, ike and the corresponding channel-tap matrix H of the 


(22) 


code matrix Gj is given by 
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ee | (23) 
h, —h, 


Assume that the channel taps 4, and hz are obtained from perfect channel 


estimation. The receiver uses the matrix H to perform MRC as follows: 


pd-[2 thy 





h 














to 





a : (TA ILs]+ IJ} 


h 














+l ly 


lA 














= fof] AL 


“Sean [ 
2} Le"82 1 2 
S, 1 | Ajn,+ | 

+—| * ; 
S, | ||| 237, —h,n, 








= [al 











(24) 


Similarly, one can perform the same operation to demodulate s,and s, by replacing s,and 
s, with s,and s,, 4; and h2 with h3 and h4, and _X; and X2 with X3 and_X4. The resulting 


matrix is 

















(25) 


s 1 | Ayn, +h,n 
ala EL 
S4 | | gl, — Ny 


The block diagrams [5] for the transmitters and receiver are shown in Figure 7 and 8, 


respectively. 
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Figure 7. Modulator circuit from the base station for the Alamouti code. 
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Figure 8. | Demodulator circuit at the subscriber station for the Alamouti code. 


There are four transmitting antennas on the BS side and one receiving antenna on 
the SS side. The only difference is that the modulator circuit can be repeated and 
reproduced. The modulator design is now simplified as the transmit diversity encoder 


only has to process two symbols instead of four. 


2. Matrix A> (Variation of Matrix A) 


Consider the STBC defined by the matrix A> in the 802.16 standard for ORDMA 
(shown in Appendix E): 


T3 


t+ tt 


S 1 _—Ss 2 0 0 ——> Antennal 


0 0 SS;  —S4 ——} Antenna 2 


S 2 Ss : 0 ) ———>>_ Antenna3 


0 0 S4 S3 ———_ Antenna 4. 


(26) 


From the matrix shown in (26), at time slot 1 (T1), symbols s and s, are sent 
through Antenna | and Antenna 3, respectively. At time slot 2 (T2), symbols —s, and s, 


are sent through Antenna | and Antenna 3, respectively. At time slot 3 (T3), symbols 


s,and s, are sent through Antenna 2 and Antenna 4, respectively, and symbols —s, and 


s, are sent through Antenna 2 and Antenna 4 at time slot 4 (T4). 


We use MRC to decode matrix A». The first step is to prove if matrix A} is 
orthogonal. Matrix 4} is orthogonal only if (4))' 4] = 4,4) =1, where A); represents 
the transpose of matrix 42, A, represents the conjugate of 42, and (4;)° represents the 


conjugate transpose of 4; , which is the same as A, 


ss -s, 0 0 
0 O s, -s 
A) A= oo Not 
(A, ) = A, Be: Ge (27) 


Therefore, 
5,8, +858, 0 0 0 
; ’ 0 S58, + S48, 0 0 
(A])' Ay = 4,4; = er ited (28) 
0 0 S15) +5555 0 
0 0 0 S58, +848, 
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and 


SS, +858, 0 0 0 
My Sate 2 0 5,5, + 555, 0 0 
comes i 0 0 5,8, + S48, 0 (29) 
0 0 0 5,8, + S48, 


We can see that Equations (28) and (29) are quasi-orthogonal. Consequently, one 
can use the same modulator and demodulator as discussed for matrix A for this matrix. 
3. Matrix A3 (Variation of Matrix A) 


Consider the STBC defined by matrix A; in the 802.16 standard for OFDMA 
(shown in Appendix E): 


Th a. Se ae 
fy Ale, IT ot 
So S5 0 O ————» Antenna 1 
A 0 Os, -S, —______» Antenna 2 
10. 0 5 8, —____» Antenna 3 
So s, 0 0 ——_ Antenna 4. 


(30) 


From matrix given by (30), we see that in time slot 1 (T1), symbols s,and s, are 
sent through Antenna | and Antenna 4, respectively. In time slot 2 (T2), symbols —s, and 
s, are sent through Antenna | and Antenna 4, respectively. In time slot 3 (T3), symbols 
s,and s, are sent through Antenna 2 and Antenna 3, respectively, and symbols —s, and 


s, are sent through Antenna 2 and Antenna 3 in time slot 4 (T4). 
We use MRC to decode matrix 43. The first step is to prove if matrix A} is 
orthogonal. Matrix A; is orthogonal only if (4))' 4} = 4,4; =/, where A] represents 


2 


the transpose of matrix 43, A,represents the conjugate of A3, and (4, )’ represents the 


conjugate transpose of 4; , which is the same as 4,: 


(A; ) = 4, = : : G1) 


Therefore, 
5,8, + 558, 0 0 0 
. 0 , + S484 0 0 
(Aj) AP = A, Ay = eee te (32) 
0 0 $45, 1 S58; 0 
0 0 0 5S, + SS, 
and 
SS, +558, 0 0 0 
¢ 0 SS, +558, 0 0 
4; (A) = 4; 4; = Tae st ae 
0 0 S85 Sg 0 (33) 
0 0 0 5,8, + S48, 


We can see that Equations (32) and (33) are quasi-orthogonal. As a result, one can 
use the same modulator and demodulator as discussed for matrix A for this matrix. Matrix 


A3, A,A} # A) A, is similar to A>, as it is quasi-orthogonal. As a result, one can use the 


same modulator and demodulator as discussed for matrix A. 


4. Matrix D (Variation of Matrix A Using Two Frequency Bands) 


Consider the STBC defined by matrix 4 and given in the 802.16 standard: 
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—> Ant] 


Frequency 
5 5 * 0 0 Ant2 Band 1 
D = 2 1 W, 
* 
0. +0. 8 <8 ———> Antl 
3 4 Frequency 
0 0 s 5. r Ant 2 Band 2 
4 2 W, 
(34) 


In this case, we utilize two frequency bands but use two fewer time slots and two 
fewer transmitting antennas. We can see from the matrix given by (34) that in time slot 1 
(T1), symbols s,and s, are sent using frequency band W, through Antenna | and Antenna 
2, respectively. In the same time slot, symbols s,and s, are sent through Antenna | and 
Antenna 2, respectively, using frequency band W. In time slot 2 (T2), symbols —s; and 
s, are sent using frequency band W, through Antenna | and Antenna 2, respectively. In 
the same time slot, symbols —s, and s, are sent through Antenna 1 and Antenna 2, 


respectively, using frequency band W 2. The above matrix D has a diversity of two, since 
the same symbol is transmitted using two antennas. For example, symbol s; is transmitted 
by Antenna | in time slot 1 (T1) and Antenna 2 in time slot 2 (T2) by the same frequency 
band, W;. Matrix D has a code rate of two, since four symbols are transmitted in two time 


slots. The transmission process is illustrated in Figure 9. 
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Figure 9. Transmission process for matrix D. 


Observe from Figure 9 that one can reduce the matrix D to 


S; Sy 

eS 

S S 

D _ 2, 1 

‘ a4 (35) 

3 4 

a 

Bio 8 


and 
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0 


* * * 
SiS) F885. F885 FSS, 


Sy> “TS 
* 
f-ae SS 
Ty es 2 1 
HS; 
* 
S, 8; 
Therefore, 
* * * * * * 
SiS, ssa 0 535, P58, Sp 838, 
; SiO B85 ey KH Soy S48 
T\* nT _ T vl 22. 372: 1" 4 “3 4° 2 
(DD SRD =). gmk a an (37) 
SiS, ESS “BGS SH SS Sa8a FS Sy 0 
* * * * * * 
SS Oey SS ESS, 0 S58, +548, 
and 
[SS # S585 F835, FS, 8 
TAHT oi ph lor eee tT Sat Saha 
D(D’) -0'D =| » |-G8) 


From (37) and (38) we see that the matrix D is quasi-orthogonal and can be 


demodulated separately using two frequency bands, W; and W2. 


In frequency band W,, 


L: ehebee 
—S, 5 | Ay ny 
fale 

—syh,+ s,h, Ny 
Mapping this to [Y] = [H][s] + [7], we get 


y= fi ee Se Val te 
h, —h, S>5 nN, 


(39) 


(40) 


The next step is to use the channel taps [/, /z] obtained via the channel estimation 


process to estimate h,; and ho, The receiver can then use [H] to perform MRC. For 
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simplicity, assume perfect channel estimation, where [h; /2] [1 ho] and apply MRC 


on Y to obtain the corresponding decision statistics vector LX], where 


od-|y Gb 


ac IIs]+IV} 
_lal +p 
lA 











ae“L yy 


lA 


= alts] Lo 


al 


“wl ialit Le. 


=|h I; i Tila in, +h, ") 
|All | 57, — A, 0, (41) 


and ||A|| is defined as 


l= lal +A)’. (42) 


nH 
_| $34, + 84h, ff 
7 -s,h,+s,h,| |n,| 


Mapping this to [Y] = [A ][s] + [7], we get 
h hy || 83 n, 
E ale | Sa 
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In frequency band W2, 


(43) 


The corresponding decision statistics vector LX] is 















































mw ye 
_jalts a ZL 
-lat) +L om 





“| ale 

S, | Al | 4; h, || n5 

= |a| 83 + ane 
S, | } ||| 437, hn, 














(45) 
From Equations (41) and (45), we see that the two noise components 7; and n2 are 
independent since E (Aim +h,n, )(Abn, —hn, ) | =(. 


The above scheme uses two less transmit antennas and two less time slots to 
transmit the same data but uses two frequency bands. This scheme is useful for high 


speed data rate communications since it can achieve a code rate of two. 
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The block diagrams [5] for the transmitters and receiver are shown in Figure 10 


and 11, respectively. 
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Figure 10. Modulator for matrix D. 
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Combiner channel likelihood 
Demod Ratios 


Diversity Sub- Log- decoder 
Combiner channel likelihood 


Demod Ratios 


Figure 11. Demodulator for matrix D. 


There are two transmitting antennas on the BS side and one receive antenna on 
the SS side. In this case, the receiver will have higher complexity compared to the 
receiver for matrix A since it receives four symbols at a time slot from two different 
frequency bands. To demodulate the received signal, two matched filters are used to 


separate the signal and process them separately, as seen in Figure 11. 


B. PERFORMANCE EVALUATION 


All of the matrices discussed in this section are quasi-orthogonal, and MRC can 
be used for detection. Therefore, one can use closed form equations to evaluate their 


performance [4]. 
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Assume a Rayleigh fading channel. The independent, identically distributed (iid) 


subchannel tap magnitudes |h;|’s are Rayleigh distributed with normalized mean square 


value E (\a,|")= 1. The combined channel tap magnitude ||’ = sla) is a y’- density 
[= 


function with 2Z degrees of freedom. The corresponding bit error probability from [2] is 








given by 
ea we a 
Px 1- 
, ee el a 
: 4 (46) 
N, [+ 4) LAT. 141 (R) 
~ 210g, M 2 = 2 
where 


aSNR 
=,|———. (47) 
1+ aSNR 


For QPSK, the minimum Euclidean distance is d,,. = ./2€, where the diversity symbol 
energy is €,. With o? = N,/2, one obtains aSNR = &,/2N,=E,/2LN,. Furthermore, 


N,=2,M =4; therefore, the bit error probability for QPSK from [4] is 


. af iS (F <4 ) 
ae a! 


For MQAM, the minimum Euclidean distance is d,,, = 6, /(M —1). Using 








(48) 


o°=N,/2, one obtains aSNR =3(e / N,)/2(M -1) =3(E, / LN,)/2(M -1). 
Furthermore, V, =4-4/ JM ; therefore, the bit error probability for MQAM from [4] is 


given by 
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pw 2r2lVM | (2k (1-0) 
: log, M ML k 4 


Ae (ee) Sei ee), 


log, M 1=0 


(49) 


The two closed forms plots for QPSK and 16QAM are illustrated in Figure 12. The 
MATLAB code is attached in Appendix A. 


Graph of Bit Error Rate vs E,/N, for QPSK 


in] —:—-- QPSK Closed Form 
—— 16QAM Closed Form 


Bit Error Rate 





E,/N, (dB) 


Figure 12. Graph of bit error rate vs E,/N,(dB) for QPSK and 16QAM. 
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C. CHAPTER SUMMARY 


The variations of the matrix A are quasi-orthogonal. Therefore, one can use a 
simple modulator and demodulator to implement the design. In addition, one can utilize 
two frequency bands to transmit at a faster rate at the expense of using two frequency 
bands. In the next chapter, we will evaluate matrix B where eight symbols are transmitted 


using four transmitting antennas and four time slots. 
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IV. IMPLEMENTATION AND PERFORMANCE EVALUATION 
OF WIMAX NON-ORTHOGONAL STC FOR OFDMA 


In this section, we shall look at matrix B defined in the 802.16 standard to 
evaluate the STC and its variations. One can assume there are multiple antennas on the 
BS side and one receiving antenna on the SS side. 

A. MATRIX B IN 802.16 

I Matrix B 


Consider the STBC defined by matrix B given in the 802.16 standard for OFDMA 
(shown in Appendix E): 


Tl T2 T3 T4 
* * 
Sy —S> Ss —S, ——___ Antennal 
B _ S5 5; Sx —S, ———>_ Antenna2 
S3 Sy S5 Ss ——»  Antenna3 


Antenna 4. 


We can see from the matrix shown in (50) that in time slot 1 (T1), symbols 


S,,8,,8,and s, are sent through Antenna 1, Antenna 2, Antenna 3 and Antenna 4, 
respectively. In time slot 2 (T2), symbols —s;,s;,—s,ands, are sent through Antenna 1, 
Antenna 2, Antenna 3 and Antenna 4, respectively. In time slot 3 (T3), symbols s,,s,,5, 
and s, are sent through Antenna 1, Antenna 2, Antenna 3 and Antenna 4, respectively, 
and in time slot 4 (T4), symbols—s; ,—s, ,s; ands, are sent through Antenna 1, Antenna 2, 


Antenna 3 and Antenna 4, respectively. This matrix uses four time slots and four 
antennas to transmit eight symbols. Use of matrix B results in a diversity of two since the 


same symbol is transmitted twice using two antennas. For example, symbol s, is 


transmitted by Antenna | in timeslot 1 (T1) and Antenna 2 in time slot 2 (T2). Use of 


3] 


matrix B also results in a code rate of two since eight symbols are transmitted in four time 


slots. The transmission process at different time slots is illustrated in Figure 13. 


Antl 


Ant2 






Time Slot 2 


Ant! 


Time Slot 4 


Time Slot 3 \/ 
RF 


Figure 13. Transmission process for matrix B. 





We use MRC to decode matrix B. The first step is to determine if matrix B is 


orthogonal. Matrix B is orthogonal if and only if (B')’B'=BB"'=I1, where 
B' represents the transpose of matrix B, B represents the conjugate of B, and (B")’ 


represents the conjugate transpose of B" , which is the same as B : 
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(B") = B = 


Therefore, 


* * 
SS; + S,5S>5 
* 


* 
+5585 + SS, 
* * 
SS e ESS 


(BY B= 2 B= 


S\S4 — $385 
* 
5.8. F885 


and 


* * 
Sy S84 Ss 
* * 


* * 
SoS. 88, 


* * 
SS, + S585 


* * 
TSS 6 + SoSg 


* * 


S55 = 5,8, 


* * 
SiS Sa 


* * 
S587 S58, 
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* * 
555) F558, 


ef 
S85 08, 


* * 
So Ses 


* * 
S55) S85 


* * 
+S.So+ SoS, 


* * 
S655 + S782 


(51) 


* 


* 
SS) 858, 


* * 
S786 + SeS5 


* * 
S)S, FSyS> 


(52) 


* * 
S555 +55; 


* * 
S385 55) 





* * 
TSS + SgSg | 


* * * * 
SS, + S585 S)S5 + S5S6 S35, — S)S> 


* * 0 * * 
+535; +848, +5357 + S48, “—S5Se PSS 
* * * * OK * * 
‘ Si Sys. Dee See Se es, 
* OR 


* * * * * 

+53,5, +5,8, +5,5,—S,S, +88, — SS; 

B'(B') =B'B = 
S58) + S6S5 S)S6 — S585 S5S'5 os S656 


* * * * 
TSS PSS “ASSESS... SSS PSS, 





* * OK * * * * 
Si ose SHS Silly 0 SSeS 8 
* OO * OK * * * * 
S885 HS 48s° ASS Sys +5757 + 5,5, | 


From (52) and (53), we can see that matrix B is non-orthogonal. We can derive 
T 
the MRC equation to demodulate Ss = [s, S, S, S, S;, Se S, S| directly 
from the four antennas. This requires the performance of MRC on matrix B with the 
T 
receive matrix [Z] =[G][h] + [n], where n=|n,, Ny aie n,, | is the Gaussian noise 


vector. The matrix B' is the G matrix in this case. Therefore, 


SS, 8, Sy |, i 
* * * * h 
gu|2 8 84 83 ||| 
Bes Ws ie ee I) | ae 
* * * * h 
Sep es: Bae Segall MN, 
(54) 
Sh, +s,h, +s;h,+8,h, n, 


—s,h,+s,h, —s,h, +s,h, ye 
S;h,+S,h, +s,h, +5s,h, n, 


* * * * 
—sh,—sh,+s;h,+sh,| | M4 
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The vector Z=Gh+n can be mapped into Y=Hs+N by changing the second 
element of Z=[Z, Z, Z, Zi , which is Z, =—s3h,+s,h, —s,h,+5,h,, into its complex 
conjugate Z;=-s,h+s,h,—s,h,+s,h,, and the fourth element, which is from 


Z, = —8,h,— s,h, + s;h,+s,h,, into its complex conjugate Z, =—s,h, —s,h,+s,h,+5,h,. 


Thus, 
LZ, Sh+ Sh, + 83h, + Sh, n, 
Y=As+n= 2, = —syh, + S| h, 7 syhy F shy n, 
ZL, S;h,+ Sh, + sh, + S,h, n, 
Z, —s,h, - sh, + sh, + Shy n, 
re 
. (55) 
h, h, h, h, 0 0 0 0 Sy n, 
h = * h -h 0 0 0 0 é ri 
= 2 h, 4 3 4 + 2 =[A][s]+[n]. 
” ? . 0 h, hy h, h, Ss n, 
: hy hy hy hy || 8, ny 
Sy 





From Equation (55), there are eight unknown symbols s,,5,,5,,5,,5,,5, and s,. 


39555 
Since the noise components are not independent, the eight equations are not enough to 
accurately decode the symbols even if one knows (1, fo, h3 and h4 from the channel 
estimation. We require at least two copies of the same data to accurately decode the 
symbols. The received Z matrix for two receiving antennas is shown in (56), and the 
corresponding Y matrix is shown in (57). The MATLAB code attached in Appendix C is 
used to simulate the results for both one receiving antenna and two receiving antennas. 


For the two receiving antenna case, 


4] 











hy 
hy 
SS, Ss Se Sy Sy Sg | Ag ai 
a Ss; SS; 8, S; 8, 8, Sy 55 || Ih se 
Ss SS, Sy Ss Sg Sy Sy || yy is 
7 ‘ * 5 7 7 . : 
—S, Sg S$; Sg 8; Sy S83 56 || 4» Nn, 56 
ies (56) 
LIo4 | 
Shy + Shy, + S3hy3 + Syhyy + Shy + SyMyy + Szhy, + 84ly, nN 
e —sh,, + s,hy, me + ssh, Se hye shy — Syhy; PSs las " Us 
Shy, + Schyy + S7hy3 + Sghyy + Sshy + Solyy + S7My3 + Shy, nN, 
8) Shh a Solis + seg — Sqhy— S.Ny, me sy, # Silyy Ny 
and 
Z, SM + Syhyy + S3hy3 + Syhyy + Shy, + SyMyy + Syhy3 + Syhy, nN 
vee eye Z, a —s,h,,+ shy, — Ss - syIi4 — syhy,+ shy, — Sify, + S3Ing j n, 
Z; S5h,,+ Sh, + Sh, + Shy, + 8;h,,+ Shy, + Sh); + Shy, Ny 
Z, —s5h,, = Sli a ssh; + Selig Fe Shy; - Sh ate sshy, ats Sel, nN, 
fe 
So (57) 
Ay hy My My hy hy hy My | 8; nh 
= hy hy, his —hy hy; hy, hy, —hy, S4 i ny 
Ay My hy My My My hy My |] Ss ny; 
hy, hy hy, -h,; bi his -h, 1 —h, S6 ny 
Sy 





The block diagrams [5] for the four transmitters and two receivers design are 


illustrated in Figure 14 and 15, respectively. 
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Figure 14. Modulator circuit from the base station for matrix B. 






Diversity Sub- Log- decoder 
Combiner channel likelihood 
Demod Ratios 


Diversity Sub- Log- decoder 
Combiner channel likelihood 


Demod Ratios 


Figure 15. Demodulator circuit at the subscriber station for matrix B with two 
receiving antennas. 


There are four transmitting antennas on the BS side and two receiving antennas 
on the SS side. The two receivers at the SS provide two copies of the message to 


accurately demodulate the received symbols. 
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2. Matrix B, (Variation of Matrix B) 


Consider the STBC defined by the matrix B, given in the 802.16 standard (shown 
in Appendix E): 


Tl T2 T3 T4 
S —s, Ss —s, ——__ Antennal 
B = S5 5, Sc —S, —__ Antenna2 
2 * * (58) 
S 4 S 3 S 8 S 6 ———» _ Antenna3 
* * 
S3 Sy Sy Ss ——> Antenna 4. 


We can see from matrix in (58) that in time slot 1 (T1), symbols s,,s,s, and s, 
are sent through Antenna 1, Antenna 2, Antenna 3 and Antenna 4, respectively. In time 
slot 2 (T2), symbols —s,,s,,—-s,ands, are sent through Antenna 1, Antenna 2, Antenna 3 
and Antenna 4, respectively. In time slot 3 (T3), symbols s,,s,,s,and s, are sent through 


Antenna 1, Antenna 2, Antenna 3 and Antenna 4, respectively, and in time slot 4 (T4), 


symbols—s_,—s,, s,ands, are sent through Antenna 1, Antenna 2, Antenna 3 and Antenna 


4, respectively. Use of this matrix requires four time slots and four antennas to transmit 
eight symbols and results in a diversity of two since the same symbol is transmitted using 


two antennas. For example, symbol s, is transmitted by Antenna | in time slot 1 (T1) and 


Antenna 2 in time slot 2 (T2). Use of matrix B, also result in has a code rate of two since 


eight symbols are transmitted in four time slots. 

We use MRC to decode matrix B. The first step is to determine if matrix Bj is 
orthogonal. Matrix BY is orthogonal only if (B!)' BJ = BB} =1, where B; represents 
the transpose of matrix Bo, B, represents the conjugate of Bo, and (B;)" represents the 


conjugate transpose of B; , which is the same as B, - 
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Therefore, 


* * 
SS; + S585 


* 


* 
+5585 + S55, 


* * 
S<S¢ FE S254 


(Bl) By =BB} = 


* * 
SS Sys, 


* * 
S654 + S58 


* * 
Sea ass 


and 


8S, —S, Ss 
* a 
So Sy; SG 
* co 


* * 
SS; +E S5S5 


* * 
TSS 6 + SeSg 


* * 
S45, F855, 


* * 
S58, — 55 


* * 
SSS thoes 
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* * 
SS; = S58, 


* * 
S586 + SgS5 


* * 
S)55 F585 


* * 
SS aS, 


* * 
TS Se + SeSQ 


* * 
S5S6 + SeS'5 


(59) 


* * 
S38; +8584 


* * 
SS, 39,9, 


* * 
HSS S287 
(60) 
* * 
S85 85 


* * 
S655 Pilg, 


* * 
ESo54 AS GS 


* * 
S)S, + S585 


* * 
+538, + 545, 


By (B;) =B,B = 


* * 
SS; + S5S>5 


* * 
+5,8, + 5,5, 


* * 
S)S5 + S556 


* * 
Sass ToS 


* 


5) S6 


* * 
— SoS5 


* Ok 


* OO 
TSG Se SS 


Shs SS, 


SSo 858, 


* * 
S55 — 284 


* * 
+8785 F SoS; 





: , : ; (61) 
S58; + S685 S\S— — SoS S5S5 + SoS 6 
* * * * 0 
TS5S3 + SoS 4 TS3Se + S48, +S5S5 + S2S¢ 
* OK * OK * * * * 
S355 — S| S, SS, — SS 0 S.S5 + S6S¢ 
* OF * OK * * * * 
—S Se + S486 +8386 — S485 T8785 + 5,55 | 
From (60) and (61), we can see that matrix Bz is non-orthogonal. We use 
simulations to evaluate performance when it is used. 
3: Matrix E (Variation of Matrix B Using Four Frequency Bands) 
Consider the STBC defined by matrix B given in the 802.16 standard: 
Frequency 
Band 1 
Wi 
Tl TT2 1 T2 Frequency 
Band 3 
Pe le 2 
* * 
Antennal = q——_ Sy —S, | Ss —S, -—3 Antenna | 
* * 
Antenna2 ¢q—_ E S, Sy S 6 —S. ———» Antenna2 
= * * 
Antenna | <——W— S3 ~S, S, Ss pet <hotenna 2 
* * 
Antenna2 <——— Sy S 3 8 S 6 ? Antenna 1. 
Frequency 
Frequency Band 4 
Band 2 W, (62) 
W 


In this case, we utilize four frequency bands but use two fewer time slots and two 


fewer transmitting antennas. We can see from the matrix given by (62) that in time slot 1 
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(T1), symbols s,and s,, sand s,, s,and s_, s,and s, are sent by Antenna | and Antenna 
2 using frequency band Wi, W2, W3 and W4, respectively. In time slot 2 (T2), symbols 
—s, and s,, —s, and s,, —s; and s., —s, and s, are sent by Antenna 1 and Antenna 2 
using frequency band Wi, W2, W3 and W4, respectively. Use of matrix F results in a 
diversity of two since the same symbol is transmitted using two antennas. For example, 
symbol s,is transmitted by Antenna | in time slot 1 (T1) and Antenna 2 in time slot 2 
(T2) over the same frequency band W;. Use of matrix E results in a code rate of four 


since eight symbols are transmitted in two time slots. The transmission process for 


different time slots is illustrated in Figure 16. 
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Figure 16. Transmission process for matrix E. 
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Observe from Figure 16 that one can reduce the matrix E to 
E 


and 


(64) 
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fee te ie PRIS re, Sls <8 
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5,8; + 558, + S55 + 5,3, 

+5,5, + 565g + S78; + SS 0 
Be) = Be = 
5,5, + 5,5, + 8,5, + 8,5, 
+5585 +868, +848, + SpSy 
From (65) and (66), we can see that matrix FE is non-orthogonal but is quasi- 


orthogonal. One can demodulate it using four frequency bands, Wi, W2, W3 and W4. 


ral 
—3, S| A, ny 
xa} 
—s,h,+s,h, Ny 
Mapping this to [Y] = [A][s] + [7], we get 


L 1 h, Sy n 
fe SEH] “ 


In frequency band W,, 


(67) 


=~ 





=> 


The next step is to use the channel taps [/), h2] obtained via the channel estimation 
process to estimate h; and fy, The receiver can then use [H] to perform MRC. For 
simplicity, we assume perfect channel estimation, where [/,, 2] —[/1, ho] and apply 


MRC on Y to obtain the decision statistics vector LX], where 
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=f fe || 
S| Al A; —h, || n 
= || Sy 1 Ha 
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and ||A|| is defined as 


lal = vial +o). 


In frequency band W2, 


| dnl le] 
Wp 85 || AG 
2 Sh, +s,h, Z n, 
—s,h,+ s,h, Ny 


Mapping this to [Y] = [H][s] + [n], we get 


ee . ca eile ay 
hy —h, |S, Ny 
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(69) 


(70) 


(71) 


(72) 


The corresponding decision statistics vector LX] is, 


|, 


LA IIs] +[N]} 


| 
r> 
LH, 





[Y] 
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| 
n> 
LH, 





a 














aR +P LAT 


lal h 














= lofts] Lp 


lal 


“A? |g 
84] Alla: A JL 

“lf ee 
S, | ||| Am, —h, n, 

















(73) 
In frequency band W3, 
Z= Ss Sy h, + ny 
-s, s,|/h,| |n, 
a ssh + ssh, " n, 
~s,h,+s-h, Ny 
Mapping this to [Y] = [7][s] + [7], we get 


h +h, || 5; n 
role aL) a 


(74) 
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The corresponding decision statistics vector LX] is, 






























































(76) 


In frequency band W4, 
oe Se Weel 
i sh, + s,h, a n, 
—s,h,+ s,h, n, | 
Mapping this to [Y] = [A][s] + [7], we get 


hh, [56] .[m 
fy SE ” 


The corresponding decision statistics vector [X] is, 


(77) 
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(79) 


The above scheme uses two less transmit antennas and two less time slots to 
transmit the eight symbols but uses four frequency bands. It has a code rate of four since 
it transmits eight symbols using two time slots. This scheme is useful for high speed data 


rate communications since it can achieve a code rate of four. 


The block diagrams for the transmitters and receiver are illustrated in Figure 17 


and 18, respectively [5]. 
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Figure 17. Modulator for matrix E. 
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Figure 18. Demodulator for matrix E. 
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There are two transmitting antennas on the BS side and one receive antenna on 
the SS side. In this case, the receiver will have an even higher complexity compare to 
matrix D as it receives eight symbols in a time slot from four different frequency bands. 
To demodulate the received signal, four matched filters are used to separate the signals 


and process them separately, as seen in Figure 18. 


B. PERFORMANCE EVALUATION 


Since the matrix B is non-orthogonal, simulation is required to evaluate its 
performance. MATLAB is used for the simulation, and one can use the ML decoding 
method to perform the simulation for both one and two receiving antennas. To perform 
the decoding, take an example from [6]. Consider a MIMO system with Mr transmit 
antennas and Mp receive antennas. Assume the information bit stream to be transmitted is 
encoded into ST codeword of dimension Mr x 7, where T is the block length, comprising 


of data symbols drawn from a unit average energy constellations. We denote the ST 
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codeword by § = [s, s, ...s,], where s[k] = E [A] ni. [k]] is the transmitted vector 


symbol over the kth symbol period. The channel remains constant over the length of ST 


codeword. The signal model is 


ee rf] +n[k], k=1,2, 41, (80) 


where E, is the total average energy available at the transmitter over a symbol period, 1 
is the Mr x Mr channel transfer function and n[k] is the noise vector. All 7 received 


vector symbols in the codeword may be stacked together so that 


E 
Y= |—-HS+N (81) 
M, 


where Y =[y, y>...y,] and N =[n, n,... n; | are matrices of size Mp x T. 


To decode, assume the receiver uses a ML detection criterion based on perfect 


channel knowledge. The estimated codeword is 


S= arg min||Y — HS|{" 


- (82) 
k|— |—~—HS|k 
peel Asta] 





T 
=arg min > 
Ss eat 





where the minimization is performed over all admissible codewords s. An error occurs 
when the receiver mistakes a transmitted codeword for another codeword from the set of 


possible codewords. 


Referring to (57), for a two receive antenna system, we see that the H matrix for 


the first two columns of matrix B is 


(83) 
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The first and second columns of the matrix B shown in Equation (50) are used to 


transmits only s,,s,,s,and s,(transmit in time slots 1 and 2 only), while the third and 
fourth columns of the matrix B are used to transmit s,,5,,s, and s,(transmit in time slots 


3 and 4 only). One can separate them and run the simulation for the first and second 
column. Since they do not transmit in the same time slot, they do not interfere with one 
another. Thus, the results will be the same. Therefore, the MATLAB simulation for the 


decoding is based on 


2 


A” 


Nn 





af : : [E s 
Minimum Distance = min ||y — 7 Hs| 7 | , (84) 
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S, 


where y is the set of possible codewords and [5, 8,  §, 3] is the received codeword. 


Equation (84) is used to compute the minimum distance between the received 
signal and the constellations points in space by comparing it to every point and 
calculating the distance between them. 

There are twelve matrix variations to perform the simulation. The twelve matrices 
can be generated from Equation (98) of the 802.16 standard shown in Appendix E. The 
twelve matrices are defined as B,.; to Bi-12, and the corresponding H matrices are defined 


as H., to Hj.12. Therefore, 


* 
S; Sy Ay My hy Ing 
* * * * * 
ee S, 8 ie 8 h,, Ay My hy 
1 > « fo44y_1 = h h h h ? (85) 
S3  —Sy 21 22 23 24 
* * * * * 
S, 83 hy, hy, hy, —h,, 
* 
5S; Sy hy, hy Ih fy 
* * * * * 
Poe Sy Sy Fit hy —-h, hy; hy, 
ia aaa le h , i ’ (86) 
S, 83 21 22 24 23 
* * * * * 
S3;  —S4 hy, —h,, hy, hy, 
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* 


22 
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> > > > 
> 


24 


— 
N 


* 


22 


* 


> > > > 
N 


22 


* 


>= > > > > > 
oo Se Ges 


ay 
N 


* 


> > 
No x N 
N Ww 


a a eee Se? ee) 


(87) 


(88) 


(89) 


(90) 


(91) 


(92) 


Ss =, h, hy hy; his 
hs S> “54 H.,= h; Ing hy hy, (93) 
S; Sy hy hy Ing hy, 
S4 S5 hy, hy, h,, hy 
Ss —S; hy, hy hy; hy 
Biio= 4 2 eo 13 12 11 14 (04) 
Sz; Sy hy hy hy hy, 
S>5 -S; hy; hy, —h,, hy, 
Sy 3; hy, yy My hy; 
B Sy -S; Te hy, hy; h, —h, Bs 
1-11 « |o44,41 = > 
4 So hy, hy, hy, hy; ) 
53 Sp hy hy, hy —hy 
and 
Ss 8, Ay, My hy hy; 
Bias 4 2 Gee 14 12 1 B (96) 
2 —S4 hy My hy h;, 
S'; s; hy, My hy, —h 


The corresponding H matrix is presented beside the matrix B and is shown from 
(85) to (96). One can use the twelve matrix variations to perform the simulation. The H 
matrix is generated using the same method illustrated in Equation (57). The first four 
columns of the H matrix in Equation (57) are used in this case as we demodulate only 


four symbols s,,s,,s,and s,. The MATLAB code is given in Appendix B. The channel 


used is a Rayleigh channel, and the noise is white Gaussian noise with 0 dB variance. 
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Graph of Bit Error Rate vs E,/N. for QPSK 
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Figure 19. Simulation results for M1.1, H1-2, Hi-3, Hi-4 and QPSK closed form solution. 


From Figure 19, the simulations for matrix H.), H)-2, Hj.3 are relatively close and 
have a 4 dB drop in performance when compared to the QPSK closed form solution 
curve. This is expected since the code rate has increase from one to two. Hj4 has 
relatively poor performance compared to the rest. We will avoid using this variation of 


the matrix for transmission. 
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Graph of Bit Error Rate vs E,/N. for QPSK 
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Figure 20. Simulation results for Ms, H1-6, Hi-7, Hi-, and QPSK closed form solution. 


In Figure 20, the simulations for H.¢ and Hj.7 are relatively close and have a 4 
dB drop in performance when compared to the QPSK closed form solution curve. H\-s 
and H;. have relatively poor performance compared to the rest. We will avoid using 


these two variations of the matrix for transmission. 
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Graph of Bit Error Rate vs E,/N for QPSK 
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Figure 21. Simulation results for M1-9, 1-10, 1-11, Hi-i2 and QPSK closed form 
solution. 


In Figure 21, the simulation results for 11-9, Hi-19 and Hj.1; are relatively close 
and have a 4 dB drop in performance when compared to the QPSK closed form solution 
curve. H\.12 has relatively poor performance compared to the rest. We will avoid using 
this variation of the matrix for transmission. 

From the simulation results, we can see that H\.4, Hi-s, Hi., and Hj-12 have poor 
performance compared to M1, Hi-2, Ai-3 Ai, Ai-7, Hi-o Hi-io, and Hi11. We conclude 
that we should use these eight matrices for transmission. 

Now we should compare the performance using only one receiving antenna. The 
Hf matrix used is from Equation (55). The MATLAB code is attached in Appendix C. The 


result is illustrated in Figure 22. 
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Graph of Bit Error Rate vs Eb/No for QPSK 
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Figure 22. Simulation results for H\.;, A\-; with only one receiving antenna and QPSK 
closed form solution. 


The poor performance when using only one receiving antenna is illustrated in 
Figure 22. There is an 11 dB drop in performance compared to the design with two 
receiving antennas. 

Next, we evaluate the performance for 16QAM. Due to the large amount of time 
required to run all twelve H matrices, only the first four H matrices were simulated. The 


MATLAB code is attached in Appendix D. The simulated results are illustrated in Figure 
23. 
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Graph of Bit Error Rate vs E,/N. for QAM 
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Figure 23. Simulation results for M11, 41-2, Hi-3, H1-4 and 1}6QAM closed form 
solution. 


The 16QAM also has a 4 dB drop in performance when compared to the closed 
form solution, illustrated in Figure 23. Similar to QPSK, the simulations results for Hy.1, 
2, H\.3 are similar, but H 4 has poorer performance compared to Mj.1, Hi-2, Hi-3. We 
can now confirm that H).4 has a poorer performance and should not be used in the 


transmission process. 
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C. CHAPTER SUMMARY 


Matrix B is non-orthogonal. Therefore, programs using MATLAB is required to 
perform simulations to evaluate the performance. It is clear that with two receiving 
antennas, one can improve the performance by almost 11 dB. The code rate is two for 
matrix B instead of one for matrix A. Therefore, there is a 4 dB drop in performance 
when compared to the closed form solution. We discovered that four matrices out of 
twelve matrices have poorer performance, and these four variations should not be used 
for transmission. In addition, four frequency bands can be used to transmit the eight 
symbols at a code rate of 4 with the performance of the closed form solution. The 
implementation of the modulator and demodulator for matrix B and its variations were 


also discussed in this chapter. 


64 


V. CONCLUSION AND FUTURE WORKS 


A. CONCLUSION 


This research was focused on STCs and evaluated the use of matrix A and B given 
in the 802.16 standard. The main point of the research was to use MRC to demodulate the 


transmitted signal. 


It has been proven that since the matrix A and its variations are quasi-orthogonal, 
the noise is independent; therefore, we can demodulate the transmitted signal using the 
MRC decision statistics. The performance of the STBC for matrix A was evaluated using 


closed form equation. 


Using an adaptive system that combines all the variations of the original matrix A 
and keeping it quasi-orthogonal, we can achieve superior performance by having the 
receiver send information to the transmitter in order for the transmitter to makes changes 
to the STC. In addition, one can utilize more frequency bands to transmit the symbols 
using half the time slots and half the number of antennas used in the standard. The code 


rate can then be improved from one to four by using four frequency bands. 


Use of matrix B, which is non-orthogonal, results in a 4 dB drop in performance 
compared to the closed form solution for QPSK and 16QAM. Also, increasing the 
number of receiving antennas enables a better performance and is necessary, as it 
provides the receiver with more information to accurately determine the transmitted 


symbol. 


The variations of the matrix B for the twelve cases determined by the 802.16 
standard were simulated. Out of the twelve matrices, four were found to have poorer 


performance than the other eight and should not used for transmission. 


B. FUTURE WORK 


Due to the large amount of time required to run the MATLAB simulations, we did 


not simulate the performance of 64QAM or all twelve cases for 16QAM. Due to time 
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constraints, only four cases for the 1}6QAM were simulated. More data points need to be 
included in the MATLAB to generate more accurate results, and 64QAM should be 


investigated in detailed. 


Future researchers can also generate more variations of the H matrix and evaluate 
the performance obtained by using these variations. In addition, the MATLAB code can 


also be modified to perform a simulation when more receiving antennas are used. 
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APPENDIX A 


SThis MAT 
SThis MAT 


LAB 
LAB 


Code plots the closed form so 
code is provided by MAJ Ng We 





EsNO_lower = 0; 
ESNO_upper 60; 
num_ESNO_points = 20; 
EsSNO_dB linspace (ESNO_ 
Eb/NO values 











lower, 











EbNO_dB_OPSK 


ESNO_dB —- 10*1log10(2 























i 
EbNO_dB_QAM16 = ESNO_dB - 10%*1log10(4); 
EbNO_dB_QAM64 = ESNO_dB - 10%*1log10(6); 
EbNO_dB = [ESNO_dB; EbNO_dB_ OPSK; 
%& calculate theoretical QAM curves 
M=16; 
L=4; 
iSNR = 3*10.% (EbNO_dB ( )/10) /L/2/ (M-1) 
mu = sqrt (iSNR. / ( nee ay 





sumterm = 0; 
for 1=0:L-1 











sumterm = sumterm + nchoosek(L-1+l, 1)* 
end 
gam16BerTheory = (4-4/sqrt (M)) /log2(M)*((1 
semilogy (EbNO_dB_QAM16, gamlé6BerTheory, 
hold on 
%& QPSK 
L=4 
iSNR = 10.%* (EbNO_dB(1,:)/10)/L/2; 
mu = sqrt (iSNR./(1+iSNR)); 
sumterm = 0; 
for 1=0:L-1 

sumterm = sumterm + nchoosek(L-1+l, 1)* 
end 
qeskBerTheory = ((1-mu)/2) .* sumterm; 
semilogy (EbNO_dB_OPSK, nes ee ree Ye 
hold on 
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lution for QPSK and QAM 





EsNO_upper, num_ 


EbNO_dB_OQAM16; 


i Gee 


% 





ESNO_points) ; multiple 





EbNO_dB_OQAM64]; 


((1+mu)/2).%1; 

—mu) /2) .* sumterm; 
Yap ) 

((1+mu) /2).*1; 


) 


THIS PAGE INTENTIONALLY LEFT BLANK 
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APPENDIX B 


oe 


This MATLAB Code computes the BER for Hl to H12 matrix 
A large portion of this MATLAB code is provided by MAJ 





ole 


lear 
close all 


Q 





gqpsk_map = [ 0.7071 + 0.7071i 
—-0.7071 + 0.7071i 
0.7071 — 0.70714 
-0.7071 -— 0.70713] 












































for QPSK 
Ng Wei Gee 


bit_diff_map = [ 0 1 2 1 2 3 2 1] 
N = linspace (1e3,2e3, 20); % number of bits or symbols 
EsNO_lower = 0; 

EsNO_upper = 60; 

num_ESNO_points = 20; 

num_points_per_iter = 4; 

num_iter=ceil (N/num_points_per_iter/2); 
constellation_map = qpsk_map; 

c_factor = sqrt(2); 

angle = 31.7 / 180 * pi(); 

% initialise variables 

EsNO_dB = linspace (EsNO_lower,EsSNO_upper,num_EsNO_points) ; % multiple 
Eb/NO values 

num_bit_errors_l = zeros(1l, num_ESNO_points); 
num_sym_errors_l = zeros(1l, num_ESNO_points); 
num_bit_errors_2 = zeros(1l, num_ESNO_points); 
num_sym_errors_2 = zeros(1l, num_ESNO_points); 
num_bit_errors_3 = zeros(1, num_ESNO_points); 
num_sym_errors_3 = zeros(1l, num_ESNO_points); 
num_bit_errors_4 = zeros(1l, num_ESNO_points); 
num_sym_errors_4 = zeros(1l, num_ESNO_points); 
num_bit_errors_5 = zeros(1l, num_ESNO_points); 
num_sym_errors_5 = zeros(1l, num_ESNO_points); 
num_bit_errors_6 = zeros(1l, num_ESNO_points); 
num_sym_errors_6 = zeros(1l, num_ESNO_points); 
num_bit_errors_7 = zeros(1l, num_ESNO_points); 
num_sym_errors_7 = zeros(1l, num_ESNO_points); 
num_bit_errors_8 = zeros(1l, num_ESNO_points); 
num_sym_errors_8 = zeros(1l, num_ESNO_points); 
num_bit_errors_9 = zeros(1, num_ESNO_points); 
num_sym_errors_9 = zeros(1l, num_ESNO_points); 
num_bit_errors_10 = zeros(1l, num_ESNO_points); 
num_sym_errors_10 = zeros(1l, num_ESNO_points); 
num_bit_errors_11l = zeros(1l, num_ESNO_points); 
num_sym_errors_11l = zeros(1l, num_ESNO_points); 
num_bit_errors_12 = zeros(1l, num_ESNO_points); 
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num_sym_errors_12 = zeros(1, num_EsNO_points) ; 


tie 
for j_iter = 1:num_EsNO_points 
for i_iter= 1:num_iter(j_iter) 





2 


& Transmitter 





















































symbol_orig = randi([0 3],4,1) ; 

s = constellation_map(symbol_origt1)j; 

n = 1/sqrt (2)* (randn(4*2,1) + lix*randn(4*2,1)); % white 
gaussian noise, OdB variance 

h1ll = 1/sqrt (2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h12 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 
channel 

h13 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h14 = 1/sqrt (2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h21 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h22 = 1/sqrt (2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h23 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h24 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 





2 


%& instantaneous noise factior 
iNoiseFactor = 2*sqrt (2) *10% (-EsNO_dB(j_iter) /20); 





% Channel and noise Noise addition for QO 

H1 = [ hill hl2 h13 hl4 
conj(h12) -conj(h11) conj(h14) -conj(h13) 
h21 h22 h23 h24 
conj(h22) -conj(h21) conj(h24) -conj(h23)]; 


H2 = [ hil hl2 h14 hi13 
conj(h1l2) -conj(h11) conj(h13) -conj(h14) 
h21 h22 h24 h23 
conj (h22) -conj(h21) conj(h23) -conj(h24)]; 


H3 = [ hill h13 hl2 hil4 
conj(h13) -conj(h11) conj(hl4) -conj(h12) 
h21 h23 h22 h24 
conj (h23) -conj(h21) conj(h24) -—conj(h22)]; 


H4 = [ hil h1l4 hl2 hi13 
conj(h13) -conj(h11) conj(hl2) -conj(h14) 
h21 h24 h22 h23 
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conj(h23) -conj(h21) conj(h22) -conj(h24) ]; 
H5 = [ hill h13 hl4 hl2 
conj(h14)  -conj(h11) conj(h13) -conj(h12) 

h21 h23 h24 h22 
conj(h24) -conj(h21) conj(h23) -conj(h22)]; 
H6 = [ hill hl4 hl3 hi2 

conj(h14)  -conj(h11) conj(h12) -conj(h13) 

h21 h24 h23 h22 

conj(h24) -conj(h21) conj(h22) -conj(h23)]; 
H7 = [ hill hl3 hl2 hi4 

conj(h12) conj(h14) -conj(h11) -conj(h13) 

h21 h23 h22 h24 

conj(h22) conj(h24) -conj(h21) -conj(h23)]; 
H8 = [ hill h1l3 hl4 hil2 

conj(h12) conj(h13) -conj(h11) -conj(h14) 

h21 h23 h24 h22 

conj(h22) conj(h23) -conj(h21) -conj(h24)]; 
H9 = [ hill hl2 hl13 hi4 

conj(h13) conj(h14) -conj(h11) -conj(h12) 

h21 h22 h23 h24 

conj(h23) conj(h24) -conj(h21) -conj(h22)]; 
H10 = hill hl4 hl13 hl2 

conj(h13) conj(h12) -conj(h11) -conj(h14) 

h21 h24 h23 h22 

conj(h23) conj(h22) -conj(h21) -conj(h24)]; 
H11l = hll hl2 hl14 hi13 

conj(h14) conj(h13) -conj(h11) -conj(h12) 

h21 h22 h24 h23 

conj(h24) conj(h23) -conj(h21) -conj(h22)]; 
H12 = hll hl4 hl2 hi13 

conj(h14) conj(h12) -conj(h11) -conj(h13) 

h21 h24 h22 h23 

conj(h24) conj(h22) -conj(h21) -conj(h23)]; 
yl = H1*s + iNoiseFactor*n(l:num_points_per_iter) 
y2 = H2*s + iNoiseFactor*n(l:num_points_per_iter) 
y3 = H3*s + iNoiseFactor*n(l:num_points_per_iter) 
y4 = H4%*s + iNoiseFactor*n(l:num_points_per_iter) 
yo = H5*s + iNoiseFactor*n(l:num_points_per_iter) 
y6é = H6*s + iNoiseFactor*n(1l:num_points_per_iter) 
y/7 = H7*s + iNoiseFactor*n(l:num_points_per_iter) 
y8 = H8*s + iNoiseFactor*n(1l:num_points_per_iter) 
y9 = H9*s + iNoiseFactor*n(1l:num_points_per_iter) 





1 


1, 


ff 


1, 


f 


f 


1, 


, 


1, 


f 





yl0 = H10*s + iNoiseFactor*n(1l:num_points_per_iter); 
yll = H11*s + iNoiseFactor*n(1l:num_points_per_iter); 
yl2 = H12*s + iNoiseFactor*n(1l:num_points_per_iter); 
min_distl=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iterl=[sl s2 s3 s4]'; 
s_iteratel= constellation_map(s_iterl1+l1); 
distl = sum(abs(yl - H1*s_iteratel).%2); 
if distl<min_distl 
min_distl=distl; 
current_best_shatl = s_iterl; 
end 
end 
end 
end 
end 
min_dist2=1le3; 
for sl = 0:3 
for s2 = 0:3 
for so. = 0:3 
for s4 = 0:3 
s_iter2=[sl s2 s3 s4]'; 
s_iterate2= constellation_map(s_iter2+1); 
dist2 = sum(abs(y2 - H2*s_iterate2) .%2); 
if dist2<min_dist2 
min_dist2=dist2; 
current_best_shat2 = s_iter2; 
end 
end 
end 
end 
end 
min_dist3=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iter3=[sl s2 s3 s4]'; 
s_iterate3= constellation_map(s_iter3+1); 
dist3 = sum(abs(y3 - H3*s_iterate3).%2); 
if dist3<min_dist3 
min_dist3=dist3; 
current_best_shat3 = s_iter3; 
end 
end 
end 
end 
end 
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min_dist4=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 073 
for s4 = 0:3 
s_iter4=[sl s2 s3 s4]'; 
s_iterate4= constellation_map(s_iter4+1); 
dist4 = sum(abs(y4 - H4*s_iterate4) .*%2); 
if dist4<min_dist4 
min_dist4=dist4; 


current_best_shat4 = s_iter4; 
end 
end 
end 
end 
end 
min_dist5=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 


for s4 = 0:3 
s_iter5=[sl s2 s3 s4]'; 
s_iterate5= constellation_map(s_iter5+1); 
dist5 = sum(abs(y5 - H5*s_iterate5).%2); 
if dist5<min_dist5 
min_dist5=dist5; 


current_best_shat5 = s_iter5; 
end 
end 
end 
end 
end 
min_dist6=le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 


for s4 = 0:3 
s_iter6=[sl s2 s3 s4]'; 
s_iterate6= constellation_map(s_iter6+1)j; 
dist6 = sum(abs(y6 - H6*s_iterate6é) .%2); 
if dist6<min_dist6 
min_dist6=dist6; 
current_best_shat6 = s_iter6; 
end 
end 
end 
end 
end 
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min_dist7=le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iter7=[sl s2 s3 s4]'; 
s_iterate7= constellation_map(s_iter7+1); 
dist7 = sum(abs(y7 - H7*s_iterate7).%2); 
if dist7<min_dist7 
min_dist7=dist7; 


current_best_shat7 = s_iter7; 
end 
end 
end 
end 
end 
min_dist8=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iter8=[sl s2 s3 s4]'; 
s_iterate8= constellation_map(s_iter8+1)j; 
dist8 = sum(abs(y8 - H8*s_iterate8) .*%2); 
if dist8<min_dist8 
min_dist8=dist8; 
current_best_shat8 = s_iter8; 
end 
end 
end 


end 
end 


min_dist9=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 023 
for s4 = 0:3 
s_iter9=[sl s2 s3 s4]'; 
s_iterate9= constellation_map(s_iter9+1); 
dist9 = sum(abs(y9 - H9*s_iterate9) .*%2); 
if dist9<min_dist9 
min_dist9=dist9; 
current_best_shat9 = s_iter9; 
end 
end 
end 
end 
end 


min_dist10=1le3; 
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for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iterl0=[sl s2 s3 s4]'; 
s_iteratel0= constellation_map(s_iter10+1); 
dist10 = sum(abs(y10 - H10*s_iteratel0) .*2); 
if dist10<min_dist10 
min_dist10=dist10; 
current_best_shat10 = s_iterl0; 
end 
end 
end 
end 
end 
min_distll=le3; 
for sl = 0:3 
tor s2 = 033 
for s3) = 0:3 


for s4 = 0:3 
s_iterll=[sl s2 s3 s4]'; 
s_iteratell= constellation_map(s_iterl1+1); 
distll = sum(abs(yll - H11*s_iteratell).%*2); 
if distll<min_dist11l 
min_distll=dist1l; 


current_best_shat1ll = s_iterll; 
end 
end 
end 
end 
end 
min_dist12=1le3; 
for sl = 0:3 
for s2 = 0:3 
for s3 = 0:3 


for s4 = 0:3 
s_iterl2=[sl s2 s3 s4]'; 
s_iteratel2= constellation_map(s_iter12+1); 
dist1l2 = sum(abs(y1l2 - H12*s_iteratel2) .*2); 
if dist12<min_dist12 
min_dist12=dist12; 
current_best_shat12 = s_iterl12; 
end 
end 
end 
end 
end 


& count bit errors 


iis 


dummy (constellation_map (symbol 
constellation_map (current_best_shat1+1) ) 
num_bit_errors_1(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


Q 
wm 


count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat2+1) ) 
num_bit_errors_2(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


° 
oO 


count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat3+1) ) 
num_bit_errors_3(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


Q 
oO 


count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat4+1) ) 
num_bit_errors_4(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


Q 
oO 


count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat5+1) ) 
num_bit_errors_5(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 
% count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat6+1) ) 
num_bit_errors_6(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


% 


count bit errors 
dummy (constellation_map (symbo] 
constellation_map (current_best_shat7+1) ) 
num_bit_errors_7(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ) ) 4 
sum (bit_diff_map (round (abs (imag (dummy) ) ) 4 


% 


count bit errors 
dummy (constellation_map (symbol 


constellation_map (current_best_shat8+1) ) 


| origtl)- 
* c factor/2; 


num_bit_errors_1(j_iter) 


t1)) + 


tl) ); 


| origtl)- 
* c factor/2; 


num_bit_errors_2(j_iter) 


t1)) + 


tl) ); 


| origt1l)- 
* c factor/2; 


num_bit_errors_3(j_iter) 


t1)) + 


tl) ); 


| origt1l)- 
* c_ factor/2; 


num_bit_errors_4(j_iter) 


t1)) + 


Hl) ); 


| origtl)- 
* ¢ factor/2; 


num_bit_errors_5(j_iter) 


t1)) + 


tl) ); 


| origt1l)- 
* c factor/2; 


num_bit_errors_6(j_iter) 


t1)) + 


tl) ); 


| origt1l)- 
* c factor/2; 


num_bit_errors_7(j_iter) 


t1)) + 


tl) ); 


L-origth)= 
* c_factor/2; 


num_bit_errors_8(j_iter) num_bit_errors_8(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1))j; 
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Q 


%® count bit errors 
































+ 


dummy = (constellation_map(symbol_origt1l)- 
constellation_map (current_best_shat9+1)) * c_factor/2; 
num_bit_errors_9(j_iter) = num_bit_errors_9(j_iter) 
sum (bit_diff_map (round (abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
% count bit errors 
dummy = (constellation_map(symbol_origt1l)- 
constellation_map(current_best_shat10+1)) * c_factor/2; 
num_bit_errors_10(j_iter) = num_bit_errors_10(j_iter) 
sum (bit_diff_map (round(abs (real (dummy)))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
% count bit errors 
dummy = (constellation_map (symbol_orig+1)- 
constellation_map (current_best_shat11+1)) * c_factor/2; 
num_bit_errors_11(j_iter) = num_bit_errors_11(j_iter) 
sum (bit_diff_map (round(abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
% count bit errors 
dummy = (constellation_map(symbol_origt1l)- 
constellation_map (current_best_shat12+1)) * c_factor/2; 
num_bit_errors_12(j_iter) = num_bit_errors_12(j_iter) 
sum (bit_diff_map (round(abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
end 
end 
N = num_iter * num_points_per_iter; 
simBer_1 = num_bit_errors_1./N/2; % simulated ber 
simBer_2 = num_bit_errors_2./N/2; % simulated ber 
simBer_3 = num_bit_errors_3./N/2; % simulated ber 
simBer_4 = num_bit_errors_4./N/2; % simulated ber 
simBer_5 = num_bit_errors_5./N/2; % simulated ber 
simBer_6 = num_bit_errors_6./N/2; % simulated ber 
simBer_7 = num_bit_errors_7./N/2; % simulated ber 
simBer_8 = num_bit_errors_8./N/2; % simulated ber 
simBer_9 = num_bit_errors_9./N/2; % simulated ber 
simBer_10 = num_bit_errors_10./N/2; % simulated ber 
simBer_11 = num_bit_errors_11./N/2; % simulated ber 
simBer_12 = num_bit_errors_12./N/2; % simulated ber 
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APPENDIX C 


This MATLAB Code computes using Hl with only 1 receiving antenna for 
SOPSK 
SA large portion of this MATLAB code is provided by MAJ Ng Wei Gee 


clear 
close all 





gqpsk_map = [ 0.7071 + 0.7071i 
-0.7071 + 0.70712 
0.7071 — 0.70714 
-0.7071 -— 0.70713] 





bit_diff_map = [ 0 1 2 1 2 3 2 1 ] 
N = linspace (le6,2e7, 20); % number of bits or symbols 
EsNO_lower = 0; 

EsNO_upper = 60; 

num_ESNO_points = 20; 

num_points_per_iter = 4; 

num_iter=ceil (N/num_points_per_iter/2); 





constellation_map = qpsk_map; 
c_factor = sqrt(2); 
angle = 31.7 / 180 * pi(); 

















& initialise variables 

EsNO_dB = linspace (ESNO_lower,ESNO_upper,num_EsNO_points) ; % multiple 
Eb/NO values 

num_bit_errors_l = zeros(1, num_ESNO_points); 

num_sym_errors_1l = zeros(1l, num_ESNO_points); 











EL 
for j_iter = 1:num_EsNO_points 
for i_iter= 1l:num_iter(j_iter) 





Q 


% Transmitter 








symbol_orig = randi([0 3],4,1) ; 

s = constellation_map (symbol_origt1); 

n = 1/sqrt (2)* (randn(4*2,1) + li*randn(4*2,1)); % white 
gaussian noise, OdB variance 

h1ll = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 

h1l2 = 1/sqrt(2)*(randn(1,1) + 1i*randn(1,1)); % Rayleigh 
channel 

h13 = 1/sqrt(2)* (randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 
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h1l4 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 
channel 
h21 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 
h22 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 
h23 = 1/sqrt (2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 
channel 
h24 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 
channel 
%& instantaneous noise factior 
iNoiseFactor = 2*sqrt (2) *10% (-EsNO_dB(j_iter)/20); 
& Channel and noise Noise addition for QO 
H1 = [ hll h1l2 hi13 hil4 
conj(h1l2) -conj(h11) conj(hl14) -conj(h13) 
000 0 
00 0 Oj; 
yl = H1*s + iNoiseFactor*n(l:num_points_per_iter); 
min_distl=1le3; 
for sl = 033 
for s2 = 0:3 
for s3 = 0:3 
for s4 = 0:3 
s_iterl=[sl s2 s3 s4]'; 
s_iteratel= constellation_map(s_iterl+1); 
distl = sum(abs(yl - H1*s_iteratel).%2); 
if distl<min_distl 
min_distl=distl; 
current_best_shatl = s_iterl; 
end 
end 
end 
end 
end 
% count bit errors 
dummy = (constellation_map(symbol_origt1l)- 
constellation_map(current_best_shat1l+1)) * c_factor/2; 
num_bit_errors_1l(j_iter) = num_bit_errors_1l(j_iter) + 
sum (bit_diff_map (round (abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
end 
end 
N = num_iter * num_points_per_iter; 
simBer_1 = num_bit_errors_1./N/2; % simulated ber 
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APPENDIX D 


sThis MATLAB Code computes using Hl to H4 with 2 receiving antenna for 
S16QAM, due to the large amount of time needed for this code to run, 
Sthe H matrix is run separately 

SA large portion of this MATLAB code is provided by MAJ Ng Wei Gee 


clear 
close all 





gqaml6_map = [0.316227766016838 + 0.3162277660168381 











0.316227766016838 + 0.948683298050514i 
0.316227766016838 — 0.316227766016838i 
0.316227766016838 — 0.948683298050514i 
0.948683298050514 + 0.316227766016838i 
0.948683298050514 + 0.948683298050514i 
0.948683298050514 — 0.316227766016838i 
0.948683298050514 —- 0.948683298050514i 
-—0.316227766016838 + 0.316227766016838i 
—0.316227766016838 + 0.948683298050514i 
-—0.316227766016838 — 0.316227766016838i 
—0.316227766016838 —- 0.948683298050514i 
—0.948683298050514 + 0.316227766016838i 
-—0.948683298050514 + 0.948683298050514i 
—0.948683298050514 —- 0.316227766016838i 
-—0.948683298050514 —- 0.948683298050514i] 

bit_diff_map = [ 0 1 2 1 2 3 2 L. «| 

N = linspace (1le4,3e5, 20); % number of bits or symbols 

m=4 

ESNO_lower = 0; 

EsNO_upper = 40; 

num_ESNO_points = 20; 

num_points_per_iter = 4; 

num_iter=ceil (N/num_points_per_iter/m) ; 





constellation_map = gqam1l6_map; 
c_factor = sqrt(10); 
angle = 0.5*atan(0.5); 




















%& initialise variables 

EsNO_dB = linspace (ESNO_lower,ESNO_upper,num_EsNO_points) ; % multiple 
Eb/NO values 

num_bit_errors = zeros(l, num_ESNO_points); 

num_sym_errors = zeros(l, num_ESNO_points); 

tic 

for j_iter = 1:num_EsNO_points 





for i_iter= 1l:num_iter(j_iter) 
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channel 


channel 
channel 
channel 
channel 


channel 


channel 


Q 


& Transmitter 














symbol_orig = randi([0 15],4,1) ; 

s = constellation_map(symbol_orig+1)j; 

n = 1/sqrt (2)* (randn(4*2,1) + li*randn(4*2,1)); % white 
gaussian noise, OdB variance 

h1ll = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 

h1l2 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 

h13 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 

h1l4 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 

h21 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 

h22 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 

h23 = 1/sqrt(2)*(randn(1,1) + li*randn(1,1)); % Rayleigh 

h24 = 1/sqrt(2)*(randn(1,1) + 1li*randn(1,1)); % Rayleigh 


channel 


AP AP AAP NP WP AP WP AAP AP AP AP WP OP AP IP OP 


ole 





























Q 


%& instantaneous noise factior 
iNoiseFactor = 2*sqrt (2) *10%* (-EsNO_dB(j_iter)/20); 





& Channel for Hl 
H=[ hill h1l2 h13 hil4 

-conj(h12) conj (hil) -conj(h14) conj (h13) 
h21 h22 h23 h24 

-conj(h22) conj(h21) -conj(h24) con j(h23)]; 


& Channel for H2 

H=[ hill hl2 hl14 hil3 
con] (hi2) =Cony (nll) con y (his) -—cony (h14) 
HZ hea h24 B23 
con] (h22) =cony (h21) con y(h23) —eony (h24) 17 


& Channel for H3 

H = [ h1l1l h13 h12 hi14 
eong(his) =-eonj (hil) conj(h1l4) =-con ] (h12) 
h21 h23 h22 h24 


conj(h23) -conj(h21) conj(h24) -conj(h22)]; 


oe 


Channel for H4 
= [ hil hia bi2 his 

eon (hls) =conj(hill) conj(hl2) -cony (h14) 
h21- h24 h22 hH23 

conj (h23) —con7 (A211) con7z (h22) -—conj (h24) ] 


o 


y = H*s + iNoiseFactor*n(l:num_points_per_iter) ; 
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r 


min_dist=1le3; 
toe sl = Os15 
foe s2. = O15 
for s3 = 0:15 
for s4 = 0:15 
s_iter=[sl s2 s3 s4]'; 
s_iterate= constellation_map(s_iter+l); 
dist = sum(abs(y - H*s_iterate) .*2); 
if dist<min_dist 
min_dist=dist; 





current_best_shat = s_iter; 
end 
end 
end 
end 

end 

% count bit errors 

dummy = (constellation_map(symbol_origt1l)- 
constellation_map(current_best_shatt+1)) * c_factor/2; 

num_bit_errors(j_iter) = num_bit_errors(j_iter) + 
sum (bit_diff_map (round (abs (real (dummy) ))+1)) + 
sum (bit_diff_map (round (abs (imag (dummy) ))+1)); 
end 
end 
N = num_iter * num_points_per_iter; 


simBer = num_bit_errors./N/m; % simulated ber 
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APPENDIX E 


PART 16: AIR INTERFACE FOR BROADBAND WARBLESS ACCESS SYSTEMS NEEE Seo 402. 16-2009 
1 
5,37, 0 Of 5,-8, 0 0 S,-3, 0 0 
A, = [= SO a}, 4</° © S,-%d 4/9 o s,-s (97) 
o o 5, -%* SS, 0 0 0 oO S, 
oo Oo 0 SS Sy; 5, 0 0 


The mapping of subscript £ to determine the Matrix A, is given by the following formula: 
£~ mod floor (Vlagica!_data_subcarrier_nwmber_for_first fome_of code ~1W2),3)+1 
where 
logical _data_swbcarrier_number_for_first_tone_of code ~ 1,2,3,,.,, Total 4 of data subcarriers. 
‘The permuted matrix-8 (over two OFDMA symbols and two subcarriers where the first two columns refer to 


the first subcarrier over two OFDMA symbols and the last two columns refer to the second subcarrier over 
two OFDMA symbols) for 4-Tx-Rate 2 is given in six permuted matrices as shown im Equation (98). 


5, 2 S_-8"3 5-4 “ S$, “8, Sy 3, 
B, — |S Se SoS og 182 Sh) ST og |S) HAS, Ss] 
S. Sy Sy Ss 5, 5%) Sy Ste! SS Eady 
S, St, 5, Se $, 8, S Ss Sy 5) Se Mo 
; (98) 
5, 82 5", 1S) =, S, -S*) 5, ~3 4 5, Ss 
B, = [S12 Se Sel g |S Sa Sy Hs] 18 My Sy Se), 
So) Re Sy Sy Sy Sy S, Se S Ss 


Sy 5%, 5, Ss Sy St) Sy Sy 5, *) & 


The mapping of subscript £ to determine the Matrix B, is given by the following formula: 

&= mod floor (legical_data_swhearrier_mumber_for_first_tone_of code — U2), 6)+ | 
where 

ogieal_data_subearrier_menber_for_first_tone_of code ~ 1,2,3,..., Total # of data subcarriers. 
8.4,.8.3,5.1 Enhanced 4-Tx Matrix A with antenna grouping 


For 4-Tx antenna BS, transmission Matrix A in $.4,8,5,5 may be employed with adaptive antenna grouping, 
which is fed back from MS. 


When MS reports 06101110 om its CQICH, then BS shall group anteema 0 and 1 foe the first subcarrier and 
antenna 2 and 3 for the second subcarrier, In matrix form, it shall be read as shown im Equation (99). 


Copynght © 2009 (BEE. AN rights reserved ion? 
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